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Resumo 

Ainda não há uma comprovação experimental que possa validar os resultados 
obtidos da termodinâmica de buracos negros. Isso porque a radiação emitida pelo 
buraco negro, prevista pela teoria, é praticamente impossível de ser detectada devido 
ao baixo valor de sua ordem de grandeza. Na busca por indícios que possam validar a 
existência dessa radiação, o estudo de modelos anãlogos a esses objetos tem crescido 
consideravelmente nas últimas décadas. Eles permitem a idealização de experimentos 
em laboratórios que seriam impossíveis de serem realizados diretamente nos sistemas 
gravitacionais. 

Um fluido em movimento pode agir sobre o som da mesma forma que os espaço- 
tempos curvos podem influenciar na trajetória da luz na relatividade geral. Com 
isso, pode-se descrever a propagação dessas ondas sonoras através de uma métrica 
efetiva, sob a qual elas seguirão geodésicas nulas. Esta dissertação faz uma revisão 
destes estudos concentrando-se em uma analogia acústica para um buraco negro de 
Schwarzschild, demonstrando suas vantagens e limitações quando aplicada para o 
estudo da teoria de Hawking. 
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Abstract 

There isn't experimental evidences that can validate tiie results of black holes 
termodynamics. This because the radiation emitted by the black hole, predicted by 
theory, it is almost Impossible to be detected due to its low value of magnitude. In 
the search for clues that could validate the existence of this radiation, the study of 
analog models to those objetcs has grown considerably in recent decades. They allow 
the idealization of experiments in labs that would be extremely complicate to be done 
directly with the gravitational systems. 

A fluid in moviment can act on the sound the same way that curved space-tltne 
can influence on light trajectoiy in the general relativity. So, one can describe 
the propagation these sound waves through an effective metric, under wich they 
will foUow null geodesics. This thesis makes a review of these studies focusing 
in an analogy to Schwarzschild black hole using an acoustic system, showing its 
advantages and limitations when applied in the study of Hawking's theoiy. 
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Capítulo 1 



Introdução 



A física de buracos negros é uma das áreas de estudos que coloca em prova a teoria 
da relatividade geral de Einstein pois abrange sistemas com campos gravitacionais tão 
intensos que permite o surgimento de fenómenos que não são comuns no nosso dia 
a dia. A primeira idéia de buraco negro foi apresentada em 1783 pelo astrónomo 
britânico John Michell que se baseou na física Newtoniana afirmando que podiam 
existir estrelas tão compactas, com grande massa, que a velocidade de escape do 
seu campo gravitacional seria superior a da luz. Naquela época não se nomeava 
tais objetos de buracos negros e este conceito nem foi muito aceito quando foi 
implementada a descrição ondulatória da luz, que entre outras coisas, a tornava 
imune aos efeitos de qualquer campo gravitacional. 

No entanto, com a chegada do século XX, e com ele a teoria da Relatividade Geral 
de Einstein, a gravidade passa a ser interpretada como uma deformação da geometria 
do espaço-tempo, associando a ele uma curvatura que influencia na trajetória de 
qualquer partícula. Apenas alguns meses depois da publicação dos trabalhos de 
Einstein, o físico alemão Schwarzschild (2]j3) obteve sua primeira solução que assumia 
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um objeto esférico massivo na origem das coordenadas. Ele demonstrou a existência 
de uma singularidade na sua métrica onde toda a matéria poderia vir a se concentrar 
em um único ponto criando assim uma região de volume nulo e densidade infinita. 
Uma outra singularidade também surge em uma determinada região, caracterizada 
pelo que ficou conhecido com o raio de Schwarzschild, que circunda este ponto 
de densidade infinita. Porém, nem mesmo Schwarzschild tomou estas situações 
extremas como sendo fisicamente possíveis de ocorrer no universo. 

Este quadro mudou em 1939 com o trabalho do físico americano Oppenheimer 
|4|[^ que desenvolveu uma possível explicação para a natureza dessas singularidades. 
Ele utilizou das mesmas razões que Chandrazekhar descreveu o colapso gravitacional 
de uma anã branca em uma estrela de nêutrons para demonstrar que uma estrela de 
aproximadamente 3 massas solares poderia colapsar sobre ela mesmo até atingir as 
singularidades de Schwarzschild. Observou-se também que, para um observador 
externo ã estrela, levaria um tempo infinito para que a superfície se reduzisse 
até o raio de Schwarzschild e, devido a isso, chamavam este colapso de estrela 
congelada, pois era como se sua superfície parasse no tempo quando atingia o raio 
de Schwarzschild. 

Esta nomenclatura foi alterada quando, em 1958, David Finkelstein (6) introduziu 
um novo sistema de coordenadas que remove a singularidade da solução no raio de 
Schwarzschild e verificando que, do ponto de vista de um observador se movendo 
juntamente com a superfície da estrela, é possível atingir e ultrapassar o raio de 
Schwarzschild em um tempo finito. A partir daí foi implementado o conceito de 
horizonte de eventos para essa superfície que anteriormente pensava-se ser singular. 
Este horizonte na verdade age como um delimitador unidirecional da estrutura causal 
do espaço-tempo, ou seja, um observador na região interna do horizonte tem acesso 
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a qualquer evento que ocorra na região exterior mas, para um observador externo ao 
horizonte, tudo o que ocorre na região interior jamais estarã acessível a ele. Ressurgia 
assim o conceito de buraco negro na física, nome este que só começou a ser utilizado 
em 1967 pelo físico John Weeler durante uma palestra. Este nome ilustra o fato de 
que estes objetos altamente densos não seriam visíveis para nós por não conseguirem 
emitir qualquer sinal luminoso devido a existência do horizonte de eventos. 

A partir daí muitas características dos buracos negros começaram a ser 
apresentadas ã comunidade científica, bem como outras soluções das equações de 
Einstein, como o buraco negro de Kerr com rotação e o buraco negro de Reissner- 
Nordstrôm com carga elétrica (veja alguns exemplos em |12| ). Na verdade, todas as 
suas propriedades podem ser descritas pela sua massa, carga e momento angular, 
já que todo o tipo de massa ou radiação que possa interagir com ele desaparece no 
horizonte de eventos se tornando inacessível para nós. 

Mas um efeito muito interessante, e intrigante até então, pode ocorrer nas 
proximidades do horizonte de eventos de um buraco negro fazendo com que ele não 
seja totalmente 'negro': ele emitiria uma certa radiação conhecida como a radiação 
Hawking |13) . Foi o físico britânico, Stephen Hawking, quem fez esta descoberta no 
início da década de setenta quando aplicou os conceitos da mecânica quântica na 
física de buracos negros. 

No mundo quântico o vãcuo não é totalmente vazio, ele contém pares de partículas 
virtuais que se manifestam como imi fóton e um anti-fóton surgindo e se aniquilando 
muito rapidamente, respeitando o princípio de incerteza de Heinseberg. Acontece 
que, se estas partículas surgem muito próximo ao horizonte de eventos de um 
buraco negro, o campo gravitacional ali pode separã-las com a captura do anti-fóton, 
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emitindo um fóton para o infinito. Dessa forma haveria uma radiação sendo emitida 
no horizonte de eventos de um buraco negro. 

Esse fenómeno implica em um gasto de energia pelo buraco negro para separar 
as partículas do vácuo, o que acaba resultando em um processo de evaporação, ou 
seja, na diminuição do seu tamanho com o passar do tempo. Mais ainda, a radiação 
Hawking tem um espectro tal como o emitido por um corpo negro, permitindo que seja 
associado ao buraco negro uma temperatura conhecida como temperatura Hawking. 

Mesmo com toda essa gama de novas informações que veio com a teoria de 
Hawking, ela nunca teve uma verificação empírica. Isso porque a radiação prevista 
por ela é extremamente fraca, sendo praticamente impossível de ser detectada com as 
técnicas usadas pela comunidade científica atualmente. Isso se torna um problema 
maior se lembrarmos que ainda não existe uma teoria quântica da gravitação. Ou 
seja, Hawking utilizou de uma física desconhecida como base de sua teoria visto que 
a transformação da partícula virtual em real no estado de vãcuo ocorreria devido a 
ação da gravidade em um nível quântico. 

Muitos trabalhos foram desenvolvidos nessa área desde então, tanto na unificação 
da relatividade geral com a mecânica quântica quanto na investigação da radiação 
Hawking numa tentativa de obter informações acerca de sua validade teórica. Um 
novo caminho foi descoberto no início da década de oitenta pelo físico canadense 
William George Unruh |18| , quando ele notou que ondas sonoras se propagando 
em um fluido em movimento se comportavam semelhante ã luz sob a influência de 
um campo gravitacional. Se tivermos então um sistema acústico em que o fluido 
ultrapasse a velocidade do som, criaremos um horizonte de eventos acústico. Isso 
permitiu a criação de modelos acústicos análogos aos buracos negros da relatividade 
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geral e a idealização de experimentos para os fenómenos previstos pela teoria de 
Hawking. 

Esta dissertação pretende abordar os conceitos bãsicos mais relevantes para se 
compreender como esta analogia é feita na mecânica dos fluidos. Essa revisão teórica 
é feita concentrando-se em um modelo acústico que reproduza a física dos buracos 
negros de SchwEirzschild. A simplicidade proporcionada por essa escolha nos permite 
ir mais fundo nas discussões quanto a natureza da radiação Hawking. 

No capítulo 2 fazemos um estudo em cima da métrica de Schwarzschild e a 
física de buracos negros que ela envolve. Nos concentramos em discutir todas suas 
propriedades, desde a caracterização do sistema de coordenadas até a interpretação 
de suas singularidades. Com isso, analisamos as órbitas para partículas e raios de 
luz, e descrevemos os sistemas de coordenadas de Finkelstein e Kruskal até obtermos 
o conceito de buraco negro atued. 

O capítulo 3 tem como objetivo obter os mesmos resiiltados da teoria de Hawking. 
É feita imia introdução matemãtica importante sobre horizontes de Killings que nos 
permite inserir os conceitos de gravidade superflcial e temperatura Hawking até obter 

a formulação de uma teoria termodinâmica de buracos negros. 

O capítulo seguinte trata da demonstração de como, em um modelo actistico, a 
propagação das ondas sonoras em um fluido em movimento pode se assemelhar à 
luz se propagando sob influência de um campo gravitacional. Deduzimos a métrica 
acústica analisando como esta analogia pode ser estabelecida e quais são seus limites, 
flnalizando com a obtenção de uma expressão para a gravidade superficial acústica, 
e consequentemente a temperatura Hawking equivalente para o sistema acústico. 
Neste quarto capítulo fica claro que nosso modelo só consegue imitar a cinemãtica da 
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métrica de Schwarzschild, sendo assim um modelo limitado quanto à reprodução de 
uma termodinâmica análoga àquela existente na física de buracos negros. 

No quinto e último capítulo, utilizamos um modelo dilatônico para reduzir o buraco 
negro de Schwarzschild 4D à um sistema 2D (tempo + espaço) e, a partir desta nova 
configuração geométrica, obtemos um análogo acústico também bidimensional. Isso 
muda consideravelmente o quadro comentado no parágrafo anterior possibiliteindo a 
inclusão da dinâmica na nossa analogia e a definição de grandezas equivalentes à 
massa e entropia de um buraco negro para o sistema acústico. 
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Capítulo 2 

A Física de Buracos Negros 



Após a formulação da Teoria da Relatividade Geral, o físico alemão, Karl 
Schwarzschild obteve a primeira solução exata da equação de Einstein para a 
gravitação. Esse trabalho, publicado em 13 de janeiro de 1916 e, embora tenha sido 
relativamente simples por sugerir uma configuração espacial com simetria esférica 
no vácuo, pode ser considerado o início de toda física de buracos negros conhecida 
hoje em dia. Iremos estudar neste capítulo essa solução de Schwarzschild (8|- |10| , 
interpretando e analisando suas principais propriedades. 



2.1 A solução de Schwarzschild 

Como dito anteriormente, a métrica de Schwarzschild é uma solução das equações 
de campo de Einstein para o vãcuo, ou seja, sem fontes. Portanto, ela descreve o 
campo gravitacional fora de um corpo massivo esférico e sem rotação. Seu elemento 
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de linha, para um grupo particular de coordenadas, é escrito como |8|, 

^(^^_ 2G^^ (^^^^2 _(^^_ 2^/^ '\r^_ ^dO^ + sin^^d^^) . (2.1) 

Aqui nós adotamos para a assinatura dos coeficientes da métrica (+,—,—,—), tal 
que ds^ > O é um intervalo do tipo tempo e ds^ < O é um intervalo do tipo espaço. 
As coordenadas = {ct, r, 9, <p) são chamadas de coordenadas de Schwarzschild e a 
métrica correspondente g^v{x) tem as seguintes propriedades: 

• A Solução é Estática - a métrica é independente do tempo t. Há um vetor de 
Killing ^ associado com esta simetria sobre deslocamentos na coordenada t que 
tem componentes 

^^ = (1,0,0,0). (2.2) 

• Esfericamente Simétrica - como já foi dito antes, a métrica de Schwarzschild 
tem a simetria de uma esfera, e sua invariância sobre rotações no eixo z é 
representado pelo vetor de Killing com componentes 

7?^^ = (0,0,0,1). (2.3) 

• Massa M - se GM/c^r for muito pequeno, o coeficiente de dr^ no elemento de 
linha |2. 1) pode ser expandido em primeira ordem, tal que 

ds^ = (^1 _ (cdí)2 - (^1 + dr' - [dO^ + Wedcf") . (2.4) 

E esta é exatamente a forma da métrica de um campo gravitacional fraco e 
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estático com um potencial gravitacional [9], dado por 



4. = -!iM. (2.5) 



Isto nos leva a identificar a constante M na métrica Schwarzschild como sendo 
a massa total da fonte de curvatura do espaço-tempo. Portanto, a geometria fora 
de uma fonte esfericamente simétrica é caracterizada por uma única grandeza 
(sua massa total) e não como esta massa está radialmente distribuída dentro 
da fonte. Por conveniência nós definimos uma massa geométrica m que tem 
dimensão de comprimento, 

m = — ^. (2.6) 

Levando em consideração os itens descritos acima, podemos reescrever a equação 
\2. 1) como 



ds^={l-^\ [cdtf ' dr' - r\de^ - senHd<t?). (2.7) 



2.1.1 Caracterização das coordenadas 

Uma métrica pode ser representada por diversos sistemas de coordenadas, porém, 
se ela possui certas simetrias então existirá coordenadas preferenciais adaptadas 
conforme estas simetrias. 

Observando as componentes diagonais da métrica de Schwarzschild concluímos 
que, para r > 2m, = ct é tipo tempo e = r, = 9 e = (p são tipo espaço. 
Como a métrica é independente de í e não há termos cruzados com dt, ela é estática 
e í é identificado como uma coordenada que representa o tempo. A coordenada r é 
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um parâmetro radial com a propriedade de que a 2 -esfera t = constante e r = constante 
tem o elemento de linha 



de onde segue que a área da superfície dessa 2-esfera é 47rr^. Finalmente, as 
coordenadas 6 e (f> são os ângulos usuais das coordenadas esféricas, que sâo 
invariantemente definidas pela simetria esférica. Portanto, as coordenadas de 
Schwarzschild são coordenadas canónicas definidas invariantemente pelas simetrias 
presentes na métrica. 

A variação radial infinitesimal da distância própria é obtida da métrica fazendo dt = 
dO = d(j) = 0. Essa distância própria não é dr, mas sim dr/{l — 2m/r)^/^. A coordenada 
r portanto não mede distância própria. De forma similar, a coordenada t também 
não é o tempo próprio. A variação infinitesimal no tempo próprio é relacionada com a 
variação infinitesimal da coordenada t tal que ds = cdr. Fazendo então dr = d9 = d(f> = O 
no elemento de linha (2. 1) , teremos 



Note que dr = dt quando r — > +00 tal que a coordenada temporal t é interpretada como 
sendo o tempo próprio quando medida a uma distância infinita da fonte do campo 
gravitacional. 

A consequência desta diferença na medida do tempo localmente e no infinito é que 
a radiação enviada de uma certa posição de raio r sofre um desvio para o vermelho 
quando é recebida por um outro observador bem distante. Como o comprimento de 



ds^ 



(2.8) 




(2.9) 
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onda da radiação é proporcional ao período da vibração, a equação (2.9) nos diz que 

A=íl-— j Aoo, (2.10) 

para a relação entre o comprimento de onda A da radiação emitida em r e o 
comprimento de onda Aoo recebido no infinito. Para o desvio para o vermelho z, 
definido por z = (Aoo — A) /A, temos que 



l + z = 




(2.11) 



Considere agora um caso onde r < 2m. Nesta região, as coordenadas r et invertem 
seus papéis. Por exemplo, para um corpo em r < 2m, com 6, cp e t constantes, temos 
dr^ = giidr'^ > O, então dr é um intervalo tipo tempo nessa região. Isto significa que 
o corpo não poderá permanecer em repouso, ele irá cair até atingir a singularidade 
r = 0. Mais ainda, como os coeficientes da métrica dependem do tempo (r nesse 
caso), ela não é mais estática nessa região. No caso da coordenada t, ela passa a 
desempenhar um papel de coordenada espacial. 



2.2 Singularidades 

De acordo com a métrica j2.1| vemos que a métrica se torna degenerada em dois 
pontos específicos: r = 2m, que torna gu divergente, e r = O onde 500 é divergente. 
O primeiro ponto é conhecido como o raio de Schwarzschild (r^) e identifica uma 
singularidade de coordenadas. Uma singularidade deste tipo pode ser removida com 
uma transformação de coordenadas adequada, sendo assim considerada como um 
resultado da matemática utilizada para descrever o espaço-tempo em questão. 
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Isso pode ser verificado através do invariante escalar do tensor de Riemann fS], 
calculado para a métrica de Schwarzschild 



RatcdR'''''=^, (2.12) 



que é finito para r = r^. Por ser um escalar, seu valor é o mesmo em qualquer 
sistema de coordenadas. Porém, esse valor diverge para r = O o que implica 
que a singularidade neste tempo não é removível como a anterior. Esta então é 
uma singularidade real e um problema físico proveniente da própria métrica de 
Schwarzschild. É importante notar também que na hiperfície r = há um desvio 
para o vermelho infinito, que pode ser visto pela equação (2.1 1) . 

Se considerarmos casos onde a fonte de gravitação tenha um raio R > rg nós não 
teríamos problema algum com essa singularidade em r = 2m. Isso porque a métrica 
de Schwarzschild só é vãlida para a região exterior a fonte do campo gravitacional e 
nesse caso este ponto singular não estaria dentro do domínio da solução. No interior 
do corpo massivo o tensor momento -energia na equação de Einstein não serã nulo, 
ou seja, uma nova solução deve ser obtida para essa região. No entanto, a principal 
questão é verificar o que acontece quando a fonte tem um raio menor que o raio 
de Schwarzschild. Nesse caso a solução continuarã sendo vãlida na região exterior 
ao corpo massivo, mas passamos a ter um espaço-tempo dividido em duas regiões 
desconexas, O < r < rs e r > rs. 
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2.3 Órbitas na métrica de Schwarschild 

Vamos agora estudar as órbitas descritas por partículas de teste na geometria de 
Schwarzschild. Considerar as simetrias esférica e temporal torna esta análise mais 
simples. E para simplificar ainda mais faremos c = 1. 



2.3.1 Constantes do movimento 

As simetrias são extremamente úteis pois a partir delas podemos obter algumas 
grandezas que permanecem constantes durante o movimento da partícula. Para o 
nosso caso particular, como a métrica é independente de t e 0, as quantidades ^ • u e 
?7 • u são conservadas, onde u é a quadri-velocidade da partícula cuja as componentes 
são = dx^/dr, sendo x'^(r) o caminho seguido por essa partícula. Os vetores de 
Killing ^ e 77 são dados por \2.2) e j2.3| . Assim, temos que 

í • u = g^.u^f^u" = qqqvP = E, 



1-^)^ = ^. (2.13. 

r J ar 



onde E é uma constante e m é dado pela equação (2.6) . Se considerarmos uma 
partícula de massa unitãria, u'^ é o seu quadri-momento. Portanto, a equação |2.13) 
representa a conservação da componente temporal do quadri-momento, isto é, a 
conservação da energia. A quantidade E então, é a energia relativística por unidade 
de massa da partícula. 
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Analogamente para a simetria na direção azimutal, teremos 



• u = giJiuri''u'' = g33U'' = -L, 



rHm^e'^ = L, (2.14) 
ar 



com L sendo uma constante também. De maneira análoga, interpretamos a 
conservação da componente (j) do quadri-momento como a conservação do momento 
angular, com L sendo o momento angular por unidade de massa da partícula. 

A conservação do momento angular implica que a órbita esta contida em um único 
plano e para simplificar as contas escolhemos uma órbita no plano equatorial, onde 
9 = 7r/2, e com isso d9 = O e sen 6 = 1. Dividindo então a equação |2.7 por dr^ teremos 



utilizando as equações (2. 13 e |2. 14 escrevemos 



A quantidade Vef{r) é conhecido como o potencial efetivo e pode fornecer 



informações interessantes sobre órbitas na métrica de Schwarzschild. A figura 2.3.1 
mostra o potencial efetivo em função de r/m e nota-se que há órbitas fechadas 
e abertas, mas nada pode ser afirmado se elas são elípticas ou hiperbólicas. A 
característica mais intrigante é que nenhum valor do momento angular pode manter 
uma órbita com energia suficiente fora da região r < 2m. Pior do que isso, para 
partículas com pouco momento angular em relação a massa central, na verdade 
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L < 2\/3m, todas as órbitas acabarão dentro da região r < 2m. Isso é bem diferente do 
que é obtido para a mecânica Newtoniana, onde qualquer momento angular previne 
uma partícula de atingir o ponto r = 0. Esta é a primeira evidência que vemos de que 
para um corpo esférico com raio menor que o raio de Schwarzschild coisas peculiares 
podem acontecer para pequenos valores de r. Na verdade, isso mostra que, nenhum 
objeto massivo é capaz de sair da região r < 2m, uma vez que tenha entrado nela. 




Figura 2.1: Potencial efetivo para diferente valores de L, plotado como l/2[V^j- — 1] em 
função de r/m. Os números indicados em cada curva são os valores de L/m. 



2.3.2 Queda livre radial 

Considere uma partícula de massa de repouso mo em queda livre radial, seu 
quadri-momento é = (mon°, mou^, 0, 0) para um movimento puramente radial. De 
(2. 13 temos que 

o dt moE 

P=moj^ = J^^y (2.17) 
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e, usando ttiq = g^uP^p'^, obtemos 



p^ = ±mo{E'-l + ^y\ (2.18) 



Já vimos anteriormente que, para um observador situado bem distante da fonte 
gravitacional, a energia total por unidade de massa de uma partícula em queda 
livre radial tem um valor constante E. Por comparação, a energia e medida por 
um observador local, mantido em uma distância r fixa, com quadri-velocidade 
(n^) = (dí/dr, 0, 0, 0), é dada por n^^p^, a projeção do quadri-momento da partícula 
ao longo de n^. Para um observador em repouso no espaço-tempo de Schwarzschild, 
da métrica temos que 

dr'^ = {l-2m/r)df, (2.19) 

então 

(u^) = ([1 - 2m/r]-i/^ O, O, 0), (2.20) 



e assim 

moE 



(2.21) 



(1 _ 2^)1/2 • 

Segue que a energia medida por um observador local estacionário aumenta com a 
redução do valor de r. A equação (2.21 mostra que as duas energias são relacionadas 



por um fator de desvio para o vermelho, o que é justamente o requerido pela 
conservação da energia. 
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2.3.3 Órbitas circulares 

Uma partícula em órbita circular tem sua coordenada r constante ao longo do 
tempo, ou seja, não há velocidade nem aceleração radial. Portanto, as relações 
dr/dr = O e d^r/dr"^ = O devem ser satisfeitas. A primeira condição, juntamente com 
a equação |2. 16) , informa que E = Ve/(r). Da equação radial |2.16) também obtemos 
que 



d^r _ _l d 
dr2 ~ ^2dr 



V,'f{r). (2.22) 



Então, utilizando a condição de que o movimento da partícula não é acelerado na 
coordenada radial, temos 

dV\ dV t 

-^=2K;^ = 0, (2.23) 
dr dr 

logo, dVef/dr = 0. Portanto, órbitas circulares somente são possíveis em pontos 
extremos do potencial efetivo. 

Usando a equação |2.16) novamente, podemos obter estes pontos extremos do 
potencial efetivo. 



7-2 1 /r4 

r = — ±-J— -12L2. (2.24) 
2m 2 V 

Desta equação vemos que há duas soluções para > 12m^, o sinal negativo 
corresponderá a um máximo desse potencial, o que é um ponto instável, e o sinal 
positivo indicará um mínimo, que será estável. Quando = 12m^ há somente uma 
órbita circular e esta será ligeiramente estável, pois d^Vej/dr"^ = O, e é a órbita 'estável' 
existente de menor raio, r = 6m. Como esta órbita é somente ligeiramente estável 
qualquer perturbação, por menor que seja, levará a partícula até r = 2m. 



17 



Capítulo 2. A Física de Buracos Negros 



2.4 Órbitas de raios de luz 

o cálculo de órbitas de raios de luz na geometria de Schwarzschild é semelhante 
ao cálculo de órbitas de partículas. As quantidades ^ • u e • u são conservadas, mas 
estamos tratando agora de fótons, partículas que têm massa de repouso nula. Isso 
vai influenciar diretamente na definição da quadri-velocidade. 

Os fótons seguem geodésicas nulas por isso a curva não pode ser parametrizada 
em fiinção do tempo próprio r. Na ausência de qualquer campo gravitacional os raios 
de luz seguem trajetórias tal que x = t e esta curva pode ser parametrizada como 

= u°A, (2.25) 

onde A é o parâmetro da curva e n° = (1, 1,0,0). Com essa parametrização temos que 
= dx"/d\, ou seja, u é o quadri-vetor tangente à curva e pode ser usado exatamente 
como no caso das partículas para definir as órbitas de raios de luz. Assim, teremos 
que, tal como foi feito em |2. 13 e |2. 14 , 



2m\ dt „ , 

- = (2.26, 



rY,-L„ ,2.27, 

onde Ef e Lf são constantes. Pode-se também normalizar A tal que u coincida com o 
quadri-momento do fóton movendo-se ao longo da geodésica nula. Assim, as equações 
acima representam a conservação da energia e do momento angular para um fóton. 
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Note que, pela própria definição, u é um vetor nulo. Dessa forma. 



u • u = O 



( 



dr 
d\ 



) 



2 



2m 



) 



r 




(2.28) 



Analogamente ao caso com partículas massivas obtemos um potencial efetivo, V/(r), 
que nos ajudará a classificar melhor as órbitas de raios de luz. Na verdade, a única 
diferença entre os dois casos é o valor do produto escalar do vetor tangente (O para 
luz e 1 para partículas) . Este nível de semelhança ocorre graças à parametrização que 
fizemos. Para diferentes parâmetros a definição de u pode mudar, mas ele continua 
sendo nulo. 

O máximo do potencial efetivo ocorre quando r = 3m independentemente do valor 
de Lf. Portanto esta órbita é circular mas, devido ao fato de ser um máximo do 
potencial, ela é instável. Não há um mínimo para Vj, logo não há órbitas circulares 
estáveis para fótons. Isso nos diz que, na teoria um corpo esférico pode prender a luz 
em uma órbita circular, mas na prática isso nunca acontece. 

Devido ã liberdade na normalização do parâmetro afim. A, se o multiplicarmos 
por uma constante qualquer, a norma de u e a equação |2.28| não mudarão. Mas os 
valores áe Ef e Lf mudarão. Portanto, diferente do caso de partículas, as propriedades 
físicas das órbitas de raios de luz podem depender somente da razão entre Ef e 
Lf. Definimos então o escalar b = Lf/Ef que contém as informações acerca do 
comportamento das órbitas. Para entender melhor o papel de b considere uma órbita 
de luz no infinito, onde o espaço-tempo é plano. Podemos introduzir as coordenadas 
cartesianas e tomar este raio de luz em questão como se propagando paralelamente 
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ao eixo X, a uma distância d dele. Para r >> 2m, a quantidade b é 



ò = 



L 



E 



r^dct>/dX _^,d4_ (2.29) 



dt/d\ dt 

Para r muito grande nós temos (p ^ d/r e dr/dt —1, resultando em 



d(j) d(j) dr d 
dt dr dt 



(2.30) 



Portanto 

h = d, (2.31) 

e a constante h representa então o parâmetro de impacto de um raio de luz. 

No ponto de máximo do potencial efetivo temos Vj{^m) = L^/27m^, de onde 
podemos escrever que 

yfma. = ^,E}. (2.32) 

Assim, a característica da órbita dos raios de luz depende de ser maior ou menor 
que 27m^. Considerando órbitas vindas do infinito, se tí^ < 27m^, a luz simplesmente 
será desviada e seguirá novemiente para o infinito. Se 6^ > 27m^ a luz circulará até ser 
capturada. 

Se considerarmos órbitas que se iniciam em pequenos raios, algo entre r = 2m 
e r = 3m, no caso onde 6^ > 27m^ a luz consegue escapar. Quando < 27m^ há 
um ponto onde ela muda de sentido e é capturada pela fonte gravltacional. Isso 
indica que, se o feixe de luz começa sua órbita com momento angular pequeno, isto é, 
suficientemente próximo da direção radial, então ela escapará do campo gravltacional. 
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2.4. 1 Propagação radial da luz 



O quadri-momento de um raio de luz se propagando radialmente será 



(p") = {dt/d\,± dr/dX,0,0). 



(2.33) 



Portanto, de (2?26) , 




(1 — 2m/r) ' 



(2.34) 



onde Eoo é a energia do fóton no infinito. Como Lf é nulo, p = E^o, de modo que 



2.5 Colapso Gravitacional 

A partir do momento em que uma estrela se forma ela atinge um estado estático 
onde a energia perdida na radiaçáo é balanceada com aquela produzida pelas reações 
nucleares em seu interior, em especial a queima de hidrogénio. Esse é o estado em 
que se encontra hoje o nosso Sol, por exemplo. Porém, isso náo dura para sempre, a 
evoluçáo atingirá um estado onde náo haverá mais reações termonucleares ocorrendo, 
chegando a estrela a um estado final de equilíbrio como as anás brancas ou estrelas 
de neutrons. No entanto, para casos onde a massa é muito grande (da ordem de 3 
ou mais massas solares), a relatividade geral prevê que a pressáo interna pode náo 
suportar a atraçáo gravitacional e a superfície dessa estrela pode colapsar. 

Vamos assumir que a pressáo de suporte foi removida completamente e que a 
superfície está em queda livre. Especificamente, estamos interessados em saber 




(2.35) 
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quanto tempo leva, do ponto de vista de um observador situado a uma distância 
r » 2m fixa, para a superfície atingir o horizonte de eventos em r = 2m. 

De acordo com a equação |2.18) temos que 



^V = ^^-l + ^. (2.36) 
ar J r 



Supondo o início do movimento da partícula mo em repouso no infinito, então E"^ = 1 
e 

— =— . (2.37) 

ar y r 

Entretanto, como estamos observando a uma distância muito grande, devemos usar 
como nosso tempo próprio a coordenada t. Então, como 



ar 1 — Ira/r 



teremos 



'^^ =b /^"^ (i ^"^A j2 39) 

dt \ r \ r J ^ 

descrevendo a queda da superfície vista de um observador distante. Integramos a 
equação acima de r' = r até r' = R, onde R é muito próximo do horizonte de eventos 
(usando r' como variável de integração). 



(2.40) 



2m ) r' — 2m 



Mas a maior contribuição ocorre próximo do horizonte, então fazendo r' = 2m + e e 
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expandindo em potências de e, 



/ \ 1/2 / r, 



2m I r' — 2m e 



e assim 

í -2mln f— - iV (2.42) 

Este é o tempo gasto para que a superfície sofra um colapso até as proximidades de 
r = 2m. Note que este tempo tende ao infinito quando r — ^ 2m. 

2.5.1 Tempo de viagem da luz 

Mas o observador só consegue enxergar este evento quando os fótons provenientes 
da superfície da estrela chegarem até ele. Por isso devemos adicionar este tempo de 
viagem da luz a equação |2.42| . 

Para a propagação radial da luz, dr^ = O, e com dé* = O e (i0 = O obtém-se 

ar 1 — 2m/r 

Integrando (2^ , 

t' = r « _2mln f— - iV (2.44) 

Jj. r — 2m \2m J 

onde usou-se outra vez a contribuição maior nas proximidades do horizonte. O 
observador vê a superfície atingindo o horizonte no tempo T = t + t' , ou seja, 



r = -4mln( — -1). (2.45) 

2m 



Isto continua divergindo quando r = 2m. Portanto, para o observador distante da 
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estrela, a queda nunca termina. Ele nunca verá a superfície (ou um objeto qualquer) 
atingir o horizonte de eventos. 

2.5.2 Observador em queda junto com a superfície 

Da equação |2.37) , em termos do tempo próprio de um observador em queda livre, 
temos 




(2.46) 



representando o tempo de queda de r = fí até r = 2m. Portanto, o observador atinge 
o horizonte de eventos e o ultrapassa em um tempo finito. Isso deixa claro que o 
sistemas de coordenadas de Schwarzschild não é apropriado para descrever todo os 
espaço-tempo ao redor de uma massa esférica. Devemos utilizar um novo sistemas 
de coordenadas onde a singularidade em r = 2m seja devidamente removida. 

2.6 Diferentes sistemas de coordenadas 

Vimos anteriormente que há um grande problema na métrica de Schwarzschild 
em relação ao ponto r = 2m. Não tem como conectarmos as duas regiões que esse 
sistema de coordenadas descreve, pois a própria métrica diverge nesse ponto. Para 
isso então, utilizamos alguns sistemas de coordenadas que possibilitam descrever de 
forma contínua todo o espaço-tempo exterior ao corpo massivo. 

2.6.1 Coordenadas de Eddington-Finkelstein 

Iniciamos procurando por uma equação nas coordenadas de Schwarzschild que 
descreva raios de luz. Isto pode ser obtido integrando dr^ = O para uma propagação 
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radial. De |2.7) , temos que 



^ = dr"^ 

(l-2m/r)2" *' 



(2.47) 



onde temos definido uma nova coordenada radial de tal forma que os raios de luz 
são descritos por (i(í it r*) = O ou í it r^, = cte. Integrando esta equação obtemos 



r + 2m log 



r — 2m 



2m 



(2.48) 



Usaremos u = í + como a nova coordenada temporal, substituindo t. Com isso, 
reescrevemos a métrica de Schwarzschild e obtemos 



dr' 



2m 



dv^ — 2dvdr — r^du?' . 



(2.49) 



Com essa nova métrica, podemos obter algumas conclusões acerca da região interior 
ao horizonte de eventos, pois agora ela não é mais singular em r = 2m. Para analisar 
essas propriedades reescrevemos a equação j2.49 



como 



2drdv 



dr^+l'^-l]dv-'+r^du^ 



(2.50) 



Para um deslocamento tipo tempo ou nulo, dr^ > O, e se r < 2m o lado direito 
da equação acima será negativo. Lembrando que, feixes de luz se propagando 
radialmente, são descritos por dv = d{t + r^,) = O, e a linha do tempo de qualquer 
objeto material deve estar contida no interior do cone de luz, o deslocamento dv será 
sempre dv > 0. Assim, para fazer o lado esquerdo da equação j2.50| ser negativo 
temos que ter dr < 0. Isto implica que a singularidade r = O é um futuro inevitável 
para qualquer deslocamento na região r < 2m. 

O horizonte de eventos separa todos os eventos em seu interior do mundo exterior a 
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ele. Um explorador que supostamente atravessasse este horizonte nunca conseguiria 
voltar de sua viagem e nem sequer poderia enviar qualquer informação sobre suas 
experiências para a região exterior. E multo pior do que Isso, ele seria Incapaz de 
evitar sua própria destruição na singularidade em r = 0. Portanto, começa a ficar 
clara e mais explícita a noção de buraco negro: corpo esférico massivo com ralo menor 
que o raio de Schwarzschild. O horizonte de eventos age então como o delimitador da 
região de onde nem a luz consegue escapar. 

As coordenadas de Eddington-Flnkelsteln que usamos foram baseadas em ralos 
nulos que entram na região r = 2m no futuro (para v = cte, t aumenta ã medida que r* 
diminui). Se considerarmos u = t — r* como coordenada temporEil, descrevemos ralos 
se afastando do horizonte (í aumenta junto com r*, mantendo u constante), que na 
métrica de Schwarzschild retornará 



Isto descreve um espaço-tempo que é o inverso de um buraco negro, onde todo futuro 
de deslocamentos nulos ou tipo tempo emergem de r < 2m para o espaço ao redor. Isto 
é chamado de buraco branco. O importante é notar que, tanto buracos negros quanto 
buracos brancos são soluções matemáticas das equações de vãcuo de Einstein, mas 
isto não significa que ambos sejam fisicamente possíveis de existir. 

2.6.2 Coordenadas de Kruskal 

As coordenadas de Eddington-Flnkelsteln descrevem um espaço-tempo contendo 
imi buraco negro ou um buraco bramco, mas nunca as duas simultaneamente. Seria 
Interessante poder olhar para essas duas regiões em uma mesma métrica. Isto é 
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feito introduzindo as coordenadas de Kruskal, que também nos fornece informações 
adicionais para entender melhor a física existente no horizonte de eventos. 



Retornando a métrica de Schwarzschild, fazendo uso da equação {2.47 , obtemos 



(2.52) 



e utilizando as duas coordenadas nulas u = t — e v = t + r^, temos 



cLt"^ = f 1 I dudv — r^dõj. 



r j 



(2.53) 



Agora definimos 



V 



-e4m j 



6 , 



(2.54) 



o que nos leva à forma 



dr 



2 32m2 



'/'^"'dUdV - r'^duj, 



(2.55) 



onde consideramos r como sendo uma função de [/ e tal que 



VU = -(^ 
\2m 



(2.56) 



obtido com uso da equação equação j2.55) é a métrica de Schwarzschild 

nas coordenadas de Kruskal {U, V, 6, (p). Inicialmente a métrica é definida para U < O 
eV>0, mas através de um extensão analítica, podemos estender para U > O e V < 0. 
Note que r = 2m corresponde a UV = O, isto é, U = O ou V = 0. 

Com isso, é conveniente traçar linhas de U e V constantes para obter o diagrama 
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de espaço-tempo de Kruskal. Há quatro regiões distintas neste diagrama, cada uma 
definida de acordo com os sinais áe U e V. 




Figura 2.2: Diagrama de Kruskal. 



As regiões I e II correspondem ao espaço-tempo descrito pelas coordenadas de 
Eddington-Finkelstein para buracos negros. As regiões I e III são referentes ao espaço- 
tempo de um buraco branco de Eddigton-Finkesltein. Mas agora notamos que surge 
uma nova parte no espaço-tempo que só conseguíamos ver com as coordenadas de 
Kruskal, a região IV. As equações matemáticas então mostram a existência de um 
novo universo com a região IV, mas isto não quer dizer que ele realmente exista 
fisicamente. Ele é geometricamente idêntico ao espaço assintoticamente plano de 
Schwarzschild, a região I. 



28 



Capítulo 3 



Termodinâmica de um Buraco 
Negro 

Na década de setenta, o físico britânico Stephen Hawking publicou seu trabalho 
sobre buracos negros que se comportam como corpos negros absorvendo e emitindo 
radiação |13| . Esta radiação emitida tem um espectro de potência de Planck que nos 
permite associar ao buraco negro uma certa temperatura. 

Basicamente o que Hawking fez foi considerar os efeitos quânticos de criação 
de partículas nas proximidades do horizonte de eventos. No vãcuo, as flutuações 
quânticas se manifestam com o surgimento de um par de fóton e anti-fóton que 
se aniquilam em um intervalo de tempo que não viole o princípio da incerteza de 
Heisenberg [14]. O que ocorre é que, se estes pares surgem muito próximo do 
horizonte, um dos fótons pode ser capturado pelo buraco negro enquanto o outro 
seria emitido para o infinito. Após essa descoberta teórica, Hawking avançou ainda 
mais estabelecendo relações entre as variações de massa e ãrea com a temperatura, 
tal como é feito na termodinâmica clãssica entre energia, temperatura e entropia, 
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originando a chamada termodinâmica de buracos negros. Nas seções seguintes 
iremos, através de conceitos basicamente matemáticos, reproduzir os resultados de 
Hawking. 



3. 1 Hipersuperfícies nulas 

Se tomarmos a superfície S = r — 2m nas coordenadas de Schwarzschild e 
calcularmos o seu vetor normal W-, teremos 



lf^ = f{x)g>^'^d,S = f{x){l-^], (3.1) 



onde f{x) é uma função arbitrária. Tomando a norma desse vetor. 



r' = g,J^i'' = fix){i-^y (3.2) 



notamos que ela é nula para r = 2m. Entáo, o vetor normal ao horizonte de eventos 
tem norma nula, o que define uma hipersuperfície nula |12|. 



As hipersuperfícies nulas têm algumas propriedades interessantes. Se um dado 
vetor tangente a ela deve ser ortogonal ao seu vetor normal, entáo, 11 = implica que 
o vetor normal é, ele mesmo, tangente a hipersuperfície nula. Assim, 



, , , diX^ , , 

l"" = -7^, (3.3) 
dX 



e as curvas x'^(A) seráo geodésicas. 

Geodésicas são curvas onde o vetor tangente a elas se desloca paralelamente a ele 
mesmo, isto é, um vetor propagado paralelamente a uma geodésica é proporcional ao 
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seu vetor tangente. Com isso, podemos escrever que 



l^{V^n=9{\W, (3.4) 



/^(a^r + r^V0 = 5(A)/^ (3.5) 

dx^' 



d /dx''\ dxP 
dx'' [lã ^ dÃ 



= 5(A)Z^, (3.6) 



resultando na equação da geodésica, 



A função de proporcionalidade ^(A) pode ser escolhida de tal forma que = 
O, sendo A um parâmetro afim. 



Note que 



f]rf.^l 

^•1 = 5M-Í^r =5om^, (3.8) 



é nulo para r = 2m. Essa é a definição de um horizonte de Killing, uma superfície 
nula que tem os vetores de Killing como normaiis. Daí surge mais uma razão para 
hipersuperfície em r = 2m ser chamada de horizonte de eventos. No horizonte de 
eventos os vetores de Killing são proporcionais aos vetores normais ã hipersuperfície, 
ou seja, 

= (3.9) 
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para alguma função h. Então, 



F(y^n = o, (3.10) 



^''V^(rr)=0, (3.11) 



— Í^V f - 



Ç^V^r-T^ Ç^rô^ln|/i| =0, (3.12) 



C'^V^Ç" = <^ (3.13) 

onde K = ln\h\ é a chamada de gravidade superficial. 
Utilizando |3.9| podemos escrever que 

Ç/.V,Çp = hl^{V,h)lp + h\{VJp). (3. 14) 

Note que, do lado direito desta equação, o primeiro termo é simétrico nos índices fi e 
p, enquanto o segundo termo tem simetria nos índices u e p. Com isso, obtemos que 
sua parte anti-simétrica é nula, ou seja. 



Ç[mV.^piI=2„. = 0- (3.15) 
Como é um vetor de líilling, ele satisfaz a equação 

Ligf,u = Vu^i, + ^i,Cu = O, (3. 16) 

onde Lç é a derivada de Lie em relação ao vetor de Killing. Inserindo este resultado 
na equação |3.15 , obtemos 
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Multiplicando ambos os lados da igualdade por V^Ç^, temos que 



(3.18) 



Mas, utilizando |3.13| , podemos escrever 



P'W=2m 



P\r=2m • 



(3.19) 



Portanto, obtemos que 



r=2m 



(3.20) 



A equação (3.20 é uma forma explícita de se obter a gravidade superficial 
do horizonte de eventos através dos vetores de Killing. Utilizando a métrica de 
Schwarzschild, onde g^^ = diag{gtt, grrAA) e ^'^ = (1,0,0,0), podemos chegar a uma 
expressão da gravidade superficial em função dos elementos da métrica: 



r=2m 



r=2m 



ff^s-rr ( -g -^gtt I + g gtt\g -^gtt 



r=2m 



(3.21) 



Como na métrica de Schwarzschild ainda temos 



(3.22) 
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obtemos 

^=\^{-9tt) , (3.23) 

uma forma direta e simples para a gravidade superficial de um buraco negro 
de Schwarzschild. Substituindo a expresão para o termo temporal da métrica 
concluímos que 

/í=^, (3.24) 
ou em unidades do sistema internacional 



c3 



(3.25) 

4GM ^ ' 



3.2 Aceleração no horizonte 

Considere um corpo de massa mo inicialmente em repouso em tq. Vamos calcular 
a aceleração desse corpo supondo-o em queda livre. Seu momento é obtido das 
equações (27T7) e (27T8) 



= ( ^3^^, ^^o{E^ - 1 + 2m/r)i/2, O, O ) , (3.26) 



onde E = a/1 — 2m/ro pois p^{ro) = 0. A velocidade desse corpo, calculada em relação 
a um observador que se move lentamente em ro é dada por 



distancia própria 1 dr 

V = ; = — — . (3.27) 

tempo próprio 1 — 2m/r dt 
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Utilizando a regra da cadeia, é possível escrever a derivada acima como 



dr dr 1 dr 

-17=^0- -t: = -õ- (3.28) 

dt dr niQ dt 



Substituindo este resultado na equação |3.27| , concluímos que 

' 2m/r — 2m/ro\ '^^'^ 



1 — 2m/rQ 



(3.29) 



A aceleração medida por este mesmo observador é a variação dessa velocidade pelo 
tempo próprio, r^, tal que dt/dr = (1 — 2m/r)^^/^. Segue que 

dv dv dr dt 
dTp dr dt dTp ' 

= -^YV7^. (3.30) 
1 — Zm/ro 

Então, a aceleração para o corpo em repouso em ro será 



a{ro) = ~{1 - 2m/ro)-'/\ (3.31) 

^0 



Agora iremos considerar uma situação hipotética para obter o limite dessa 
aceleração quando r — ;> 2m. Então, imagine que um observador distante esteja 
segurando uma partícula de massa unitária por um fio ideal. A energia gasta por 
ele para suspender essa partícula de um raio r por uma distância dl é dEoo = goodl, 
onde é a métrica calculada na posição do observador. A energia da partícula 
aumenta de dEr = Qrdl, onde gr é a métrica calculada no ponto onde se encontra a 
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partícula. A conservação da energia implica que 



dEr= dE^ ^ 5oo = í?rfl-— j ; (3.32) 



mas, gr é justamente a aceleração sentida por um observador em r que foi obtida 
na equação |3.31 . Portanto, tomando o limite quando a partícula se aproxima do 



horizonte de eventos, temos que 



9oo = y*- = /í- (3.33) 

4m 



Assim sendo, a força necessãria para manter uma partícula de massa unitãria sobre 
o horizonte de eventos de um buraco negro é justamente igual ã gravidade superficial 
do buraco negro. Isto equivale a dizer que a gravidade superficial é a aceleração de 
uma partícula estática próxima do horizonte medida por um observador no infinito. 



3.3 Temperatura Hawking 

Depois dos primeiros cãlculos realizados por Hawking vãrias outras formas foram 
desenvolvidas para se obter a radiação emitida por um buraco negro (veja por exemplo 
|15| ) . Usaremos um método para investigar os efeitos quânticos nas proximidades do 
horizonte de eventos conhecido como euclideanização da métrica de Schwarzschild. 
Ele consiste em uma rotação na coordenada temporal definindo uma coordenada de 
tempo imaginãria, 

t = ít. (3.34) 
Aqui, r representa o "tempo imaginãrio" e não o tempo próprio como de costume. 
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Esta transformação nos leva à métrica de Schwarzschild euclidiana, 



que também é singular em r = 2m. Com uma transformação na coordenada radial, 



r-2m= — , (3.36) 
8m 



obtemos que, quando r — > 2m [x ^ 0), 



2m _ (kx)^ ^ 2 
r l + (/íx)2 ^ ^ ' 



(3.37) 

KX)" 

dr^ = {Kxfdx^. (3.38) 



Deste modo, teremos 



dsl « (KxfdT^ + dx'^ + -^dn^. (3.39) 



Este é o chamado espaço-tempo euclidiano de Rindler, pois é obtido através da 
euclidianização da métrica de Schwarzschild seguido da transformação de Rindler 
dada pela equação {3.36 . Este sistema de referência é centrado no horizonte de 



eventos e é simétrico sobre esta origem. 

Quando reescrevemos a equação j3.35 , tal que 



dsl = dx'^ + x^d{KT)'^, (3.40) 



ela fica idêntica ao espaço euclidiano para coordenadas polares. Porém, para 
evitarmos a singularidade cónica em x = O precisamos impor vima periodicidade sobre 
T de 2tt/k tal como é para o ângulo polar no espaço euclidiano. 
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Agora note que, no espaço euclidiano com o tempo imaginário, a integral funcional 
de um sistema quântico governado por uma hamiltoniana H é 

Z = j [D(t)\e-^^^^\ (3.41) 

onde £'[(/)] é a ação e os campos escalares 4>{x,t) são periódicos na coordenada t com 
período [ih. Com essas condições a integral funcional pode ser escrita como 

Z = tr e-^", (3.42) 

que é exatamente a função partição do sistema com uma temperatura T dada por 
/3 = l/ksT, onde ks é a. constante de Boltzman |16|. 



Aplicando estes resultados da Teoria Quântica de Campos no caso do espaço 
euclidiano de Rindler, o que equivale a escrever a periodicidade deln/K na coordenada 
r como sendo /3?i, obtemos 



ou, para h = ks = l. 



Th = ^. (3.44) 



Esta é a chamada temperatura Hawking, a temperatura de equilíbrio do horizonte 
de eventos de um buraco negro. 

Da lei de Stephan-Boltzmann para a radiação de um corpo negro temos que 



^ - -aATf,, (3.45) 



onde A é a. área do corpo negro, no nosso caso a área do horizonte de eventos. A 
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constante a é dada por 



Como 



1 TT^t"' 



E = Mc^, A = leW, (3.47) 



utilizando o valor da gravidade superficial no sistema internacional de unidades, 
temos que 



Portanto, 

dM hé 



dt G2^2 ■ 



(3.49) 



e, integrando de um tempo inicial onde o buraco negro teria massa M até um tempo 
final onde toda sua massa tenha evaporado (M = 0), isso resulta em um tempo de 
vida 

g) M'. (3.501 

Como era esperado, quanto mais massivo for o buraco negro mais tempo ele vive 
e menor é o valor da sua temperatura Hawking. Se substituirmos os valores das 
constantes em |3.43| , temos que 



Tn^j^- 1023. (3.51) 



Se um buraco negro tiver sido formado por um colapso gravitacional de uma estrela, 
sua massa seria da ordem de 10'^° Kg (ordem da massa solar por exemplo). Isso 
equivaleria a uma temperatura Hawking de 10^ K, o que é um valor muito baixo 
para ser detectado. Portanto, se existir hoje no universo buracos negros de natureza 
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estelar, eles não seriam úteis para detectar a radiação Hawking. 

Poderiamos então buscar casos onde a evaporação estaria quase no fim, com a 
massa do buraco negro sendo muito pequena. Isso possivelmente causaria uma 
explosão de radiação extremamente intensa. Mas como o tempo necessário para que 
toda a massa se evapore é proporcional a M'^, novamente vemos que buracos negros 
estelares seriam inapropriados pois seu tempo médio de vida seria da ordem de 10^^ 
anos, o que é muito maior que a própria idade do universo. 

No entanto, buracos negros com massas bem menores (da ordem de 10^^ Kg), 
poderiam se formar a peirtir da extrema densidade de matéria presente na expansão 
inicial do universo. Conhecidos como buracos negros primordiais, estes sim estariam 
hoje em seus estágios finais de evaporação. O problema aqui decorre do fato de se 
encontrar realmente estes buracos negros primordiais, visto que seu tamanho seria 
muito pequeno e com pouca influência gravitacional. 

Devemos ressaltar aqui que os cálculos originais de Hawking não levam em 
consideração a influência que a propagação dos fótons têm sobre as equações da 
mecânica quântica. Isto é válido enquanto a massa do buraco negro é maior do que 
a taxa de energia emitida pela radiação. Porém, nos estágios finais dessa evaporação, 
quando a massa é da ordem da energia emitida, seria necessário uma teoria quântica 
da gravitação consistente para realizar tais cálculos. 
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3.4 Teorema da Area 

Para um buraco negro de Schwarzschild, a área do horizonte de eventos é 
calculada a partir do elemento de área de uma esfera de raio rt'- 

dA = S-Krhdrh, (3.52) 

dA = 327:mdm, (3.53) 
G 

d A = Sn—dM, (3.54) 

CK 



OU ainda, 



dM = -^dA. (3.55) 



Substituindo a gravidade superficial pela temperatura Hawking, 



dM = TH'^dA, (3.56) 



podemos definir 



S=^A (3.57) 



e, utilizando E = Mc^, obtemos 



dE = TndS. (3.58) 



Portanto, a equação (3.56 representa a primeira lei da termodinâmica para um 



buraco negro de Schwarzschild, com a entropia dada pela equação j3.57 



Note que, enquanto o buraco negro emite essa radiação sua massa vai diminuindo, 
fato que ficou explícito quando obtemos uma taxa negativa para a variação da massa 
com o tempo. Isso diretamente implica na redução da área do horizonte de eventos ã 
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medida que o buraco negro vai evaporando e no aumento da taxa de energia irradiada 
pelo buraco. Até aí os resultados da teoria de Hawking são satisfatórios. Porém, 
quando se faz a analogia da mecânica do buraco negro com as leis da termodinâmica, 
impomos que a entropia do sistema é diretamente proporcional a área do horizonte 
de eventos. Isso implica em vima quebra imediata da segunda lei da termodinâmica, 
que diz que a variação da entropia de um sistema nunca deve ser menor que zero. 

Este impasse na formulação de uma teoria termodinâmica para buracos negros 
foi resolvido também por Hawking |17) quando ele notou que a entropia do buraco 
negro diminui com sua evaporação mas a entropia do universo externo ao horizonte 
de eventos aumenta com a radiação que chega até ele, e esse acréscimo é maior que 
a variação da entropia do buraco negro. Assim, a segunda lei da termodinâmica é 
reformulada tomando 

S = Sbn + Sext, (3.59) 

onde Sbn é a entropia do buraco negro dado por j3.57| e Sext é a entropia do universo 
exterior ao buraco negro. Com isso, a variação da entropia será sempre maior que 
zero, respeitando a segunda lei da termodinâmica. 
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Análogo Acústico a um Buraco 
Negro 

Os sistemas análogos para a relatividade geral começaram a ser desenvolvidos no 
início da década de oitenta com o físico canadense Willian G. Unruh |18| . Unruh 
demonstrou que ondas sonoras se propagando em um fluido em movimento é 
fenómeno semelhante ao da luz quando viaja em um espaço-tempo curvo. As ondas 
sonoras também sofrerão uma variação de frequência ã medida que se propagam, 
contra ou a favor, da direção de escoamento do fluido. Em especial, se o fluido atinge 
velocidades supersônicas, em um determinado ponto se formará uma barreira sonora 
tal como o horizonte de eventos de um buraco negro gravitacional. Assim, as ondas 
sonoras que estiverem dentro da região onde o fluido se movimenta mais rápido que 
a velocidade do som, nunca conseguirão escapar, definindo assim o buraco negro 
acústico. 

Mais interessante ainda é que a propagação de ondas sonoras na vizinhança desse 
buraco negro acústico pode ser descrita pelo ferramental geométrico da relatividade 
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geral através de uma métrica efetiva. Com isso, este modelo análogo aos buracos 
negros gravitacionais, se torna extremamente útil no estudo dos fenómenos existentes 
na física de buracos negros. Entre eles a radiação Hawking, principal alvo de estudo 
dessa dissertação. 

Neste capítulo iremos demonstrar como essa analogia surge nas equações 
matemáticas da mecânica dos fluidos e analisar o que ela nos revela de novo no campo 
da física de buracos negros e, consequentemente, no estudo do efeito Hawking. 

4. 1 A métrica acústica 

Com o intúito de obter uma análise mais completa possível, iremos nos concentrar 
em um modelo acústico análogo ao buraco negro gravitacional de Schwarzschild. Para 
este fim, devem ser impostas algumas considerações iniciais sobre o nosso sistema 
acústico. Trabalharemos com um fluido em movimento, tal que, em um certo ponto, 
ele atinja e supere a velocidade do som. Não iremos nos deter inicialmente a qualquer 
configuração especial que realize este feito, basta que ele ocorra. Sob estas condições, 
o fluido deve ser irrotacional para que a não haja vórtices e assim a velocidade do 
fluido possa ser descrita inteiramente por um campo escalar. Ele deve ser barotrópico 
também, ou seja, sua pressão será função somente de sua própria densidade. Isso 
garante que um fluido inicialmente irrotacional permaneça irrotacional. Por fim, para 
simplificar, consideramos o fluido não viscoso. O efeito causado quando se leva em 
consideração a viscosidade do fluido pode ser estudado em |21) . 

A dinâmica do fluido é descrito por três equações |19| : 
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A equação da continuidade 



dp 



+ V • (pv) = 0; 



(4.1) 



a equação de Euler 



(9v 



+ (v- V)v 



(4.2) 



e alguma equação barotrópica 



P = PÍP), 



(4.3) 



com p sendo a densidade do fluido, v a velocidade, p a pressão e <í> representa tanto o 
potencial gravitacional newtoniano quanto um potencial externo qualquer. 

A propagação do som é descrita matematicamente como uma perturbação linear 
das variáveis dinâmicas do meio. Desta forma, como na radiação Hawking, que 
é baseada na separação dos pares virtuais de fótons provenientes das flutuações 
quânticas do estado de vácuo, espera-se que as perturbações do fluido nas 
proximidades do horizonte de eventos acústico possa originar uma emissão de ondas 
sonoras tal como o efeito Hawking |20). 



Então, vamos linearizar estas equações com perturbações em primeira ordem em 
torno de alguma solução exata {po,Po,ipo) das equações |4.1 , |4.2) e |4.3) . Como 



escolhemos que a velocidade do fluido seja irrotacional, ou seja, V x v = O, podemos 
definir esta velocidade como 

V = - W, (4.4) 
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onde V é um campo escalar chamado potencial velocidade. Usando este potencial 
velocidade, as quantidades perturbadas serão tais que, 



p = PO + epi + 0{e^)] 
V; = ^0 + eV'i + O(e^), 



(4.5) 



onde os 0{e^) representam os termos de segunda ordem em diante, que serão 
negligenciados nos nossos cálculos. Note que estamos mantendo o potencial $ fixo, 
isto é, estamos impondo que não haja reações de fundo. Estas reações seriam a 
deformação do potencial $ causada pela propagação das ondas ipi . 



Reescrevemos as equações |4.1| - 14.3) utilizando as perturbações j4.5) , obtendo 
assim. 



dp 
~dt 



- + V • (piV^o + PoVipi) = O, 



Po 



dt 



+ VV^o • VV-i =pi, 



Pi = c,pi. 



(4.6) 
(4.7) 
(4.8) 



Para a liltima equação acima, utilizamos a definição padrão da velocidade do som, 
Cs = dp/dp, e a tomamos como uma constante. A combinação dessas três equações 
resulta em uma equação diferencial de segunda ordem para V^i 



d_ 
dt 



Po ( dil^ 



-^+vo-V^iJ 



+ V 



poVil^i 



Po 



rvo 



dt 



-+vo-VV'ij =0. (4.9) 



Esta é a equação de onda que descreve a propagação do potencial velocidade. 
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Uma vez que tenhamos determinado o valor de í/ji podemos utilizar as equações |4.6) 
e |4.7| para determinarmos pi e pi. Portanto, a equação de onda acima descreve 
completamente a propagação das perturbações acústicas. 

Agora, nosso objetivo é descrever a propagação dessas perturbações acústicas tal 
como o formalismo da relatividade geral. Para tal, introduzimos uma matriz 4x4 |21J, 
definida como 



Po 



-1 



(4.10) 



Assim, podemos reescrever a equação (4.9) como 



(4.11) 



mas o d'Alembertiano no espaço-tempo curvo |10| é definido como 



AtP = -Lô^ (V^^^^^^a.^) , (4. 12) 



onde g^^'^ é a inversa da métrica do espaço-tempo que se está calculando o 
d'Alembertiano e g é o seu determinante. Portanto, impondo a condição f^'^ = 
y/^g^'^ , obtemos nova matriz que representa uma métrica para o nosso sistema. 
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que chamaremos de métrica aciistica, 



acústica Po 



(4.13) 



Com a métrica em mãos, podemos escrever o elemento de linha acústico 



ds 



2 _ PO 



cldr + 5ij{dx' - vldt){dx^ - v^dt) 



(4.14) 



Observe que a assinatura da métrica acústica, (—,+,+,+), é exatamente como 
deve ser em uma geometria Lorentziana. É muito importante entender o que a métrica 
acústica representa no nosso sistema. Na verdade, o nosso sistema físico está situado 
em uma variedade descrita pela métrica usual de Minkowski, com a topologia básica 
do espaço K'*. Contudo, as perturbações acústicas (o campo escalar ipi) se comportam 
como se "vivessem" em uma variedade cuja métrica é dada pela equaçáo (4. 14 . 
Tanto que essa métrica acústica herda a topologia e a geometria do espaço-tempo 
de Minkowski. Devemos enfatizar também que a métrica na relatividade geral está 
relacionada a distribuiçáo de matéria no espaço-tempo, enquanto a métrica acústica 
está relacionada a distribuiçáo de matéria de uma forma algébrica |21[|22| . Mais 
corretamente, ela descreve uma espécie de distribuiçáo de velocidade do fluido, que 
age sobre o som tal como a massa sobre a luz. 
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4.1.1 Ergo-regiões e horizonte de eventos acústico 

Considere os vetores de KlUing = (1, 0, 0, 0), tal que 

9,..ee = 9u = -{(?-v''). (4.15) 

Esta expressão será positiva para ||v|| > 0. Portanto, qualquer região supersônica será 
chamada de ergo-região, o equivalente a ergo-esfera existente na geometria de Kerr 
para buracos negros com rotação. Na relatividade geral essa região é onde o espaço "se 
move" mais rápido que a velocidade da luz, sendo impossível permanecer em repouso 
lá. Na analogia acústica, tenteir ficar parado em uma ergo-região implicaria em uma 
turbulência sónica. 

Na métrica acústica podemos ver explicitamente a existência de um horizonte de 
eventos, ocorrendo quando a velocidade do fluido se iguala à velocidade local do som 
(fo = Cs). Um horizonte de eventos acústico é definido pela superfície que delimita a 
região da qual geodésicas nulas (os fônons) não podem escapar. Em toda geometria 
estacionária, como a de Schwarzschild, o horizonte de eventos é a borda da ergo- 
região. No nosso sistema acústico ocorre da mesma forma pois o horizonte separa as 
regiões subsônica e supersônica. 

Uma dos modelos mais simples para estudar os horizontes de eventos acústicos é 
de um fluido em movimento onde suas velocidades transversais são muito menores 
que a velocidade no eixo x, podendo assim ser negligenciadas. Com isso podemos 
obter o elemento de linha acústico para um fluido em movimento unidimensional 
como 

ds'^ = ^ [-{c] - vl)dt^ - 2vodtdx + dx'^ + V + dz^] , (4. 16) 

Cs 
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onde vq é o módulo da velocidade do fluido, que é um vetor no eixo x. 



4.2 Analogia à métrica de Schwarzschild 

Obtida a métrica acústica, o próximo passo é verificar se realmente podemos 
imitar a métrica de Schwarzschild e de que maneira essa analogia deve ser feita. 
Para tal, começamos com uma transformação na coordenadas de Schwarzschild para 
obter uma representação que se adapta mais a equação |4. 14) . Partimos da seguinte 
transformação na coordenada temporal [23J , 



t = ts + 



/— — _ , / V2GMr + 2Gm\ 

2V2GMr - 2GMln , 

\V2GMr - 2GM J 



(4.17) 



onde ts representa a coordenada temporal de Schwarzschild. Dessa transformação 
obtemos 

Substituindo na equação j2.7) , com c = 1, temos que 



n , 2GM\ n 2GM o o o 

ds^ = -{l dt^ + 2a/ dtdr + dr^ + r^dfí^, (4.19) 



onde é a parte angular do elemento de linha em coordenadas esféricas. Esta 
nova representação da métrica de Schwarzschild é conhecida como as coordenadas 
de Painlevè-Gullstrand, |24| , e ela cobre toda região assintótica, o horizonte futuro e 
a singularidade do buraco negro. Essa forma também é não singular no horizonte 
e transfere toda curvatura do espaço-tempo de Schwarzschild para a coordenada 
temporal, ou seja, para t = cte a parte espacial é plana. 
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Comparando o elemento de linha (4.19) com o elemento de linha acústico |4.16) 
notamos que, a menos de um fator conforme, eles são idênticos (visto que o elemento 
de linha acústico também pode ser escrito em coordenadas esféricas). Então, a 
analogia aconteceria de forma exata se escolhêssemos um fluido cuja densidade fosse 
independente da posição, tomãssemos a velocidade do som também como constante 
e definíssemos o perfil da velocidade do fluido como sendo 

vo = y/2GM/r. (4.20) 



Porém, esta configuração do nosso sistema acústico não estaria de acordo com a 
equação da continuidade (4.1 que deve ser satisfeita, pois a usamos para deduzir a 



própria métrica acústica. Isso nos leva a concluir que uma métrica acústica conforme 
com a métrica do sistema gravitacional é o mais próximo que podemos chegar com 
essa analogia. 

Contudo, isso é extremamente satisfatório para análise da radiação Hawking, 
visto que as grandezas envolvidas nessa teoria (gravidade superficial em especial) são 
invariantes perante transformações conformes da métrica (25). Este modelo então 
nos fornece um sistema munido de um horizonte de eventos que, assim como na 
relatividade geral, emite fônons segundo o argumento de Hawking. 



4.3 Gravidade superficial e Temperatura Hawking 

A radiação Hawking é caracterizada exclusivamente pela sua temperatura que por 
sua vez é obtida diretamente da gravidade superficial do buraco negro. Portanto, 
se definirmos uma gravidade superficial para o sistema acústico, conseguimos obter 
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uma equação para a temperatura Hawklng equivalente. Para isso, iremos inicialmente 
transformar a métrica acústica |4.16| para a forma tradicional de Schwarzschild 
usando para a nova coordenada temporal 



T = t + 



2vx 



2 

V X 



In 



Cs+V 



(4.21) 



tal que 



dr = dt + 



zdx. 



(4.22) 



Esta nova coordenada temporal, r, é sugerida a partir da transformação |4.17 usando 



V cx y^2GM/r e fazendo os ajustes de dimensão com c^. Substituindo 4.22 
elemento de linha acústico obtemos 



no 



ds' 



Po 

Cs 



v^)dT'^ + dx^ + dy^ + dz^ 



(4.23) 



Sob este novo sistemas de coordenadas podemos ver claramente que a geometria 
acústica é estãtica, isso devido ã ausência dos termos cruzados de tempo e 
espaço. Então, os vetores de Killing são normais ao horizonte de eventos, tal 
como foi demonstrado no capítulo anterior. Podemos utilizar a equação (3.21 para 
calcular a gravidade superficial acústica, substituindo as componentes da métrica de 
Schwarzschild pelas componentes da métrica acústica (4.23) , obtendo que 



Po 



(4.24) 



Assim, 



J_d^ 
2 Po dx 



-gttj 



2po dx 



vl) 



(4.25) 
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Note que a gravidade superficial deve ser calculada quando vq = Cs que é a definição 
do horizonte de eventos acústico. Usando a propriedade da derivada do produto de 
funções, escrevemos 



d 



-{cl 



vl) 



Cs dx 
c= dx 



dx Co 



Po d 2 2\ 

[Cs - Vq) 



VO=Cs 



Por fim, substituindo (4.26) em (4.25) , obtemos 



VO=Cs 



(4.26) 



1 d ^ 
2co dx 



(4.27) 



Com isso, utilizando a equação (3.43 obtemos a temperatura Hawking para o análogo 
acústico, 

(4.28) 



rr ^ d , 2 2\ 
T= -. : -(Cs-Vq) 



AncskB dx 



Dessa forma, o horizonte acústico emite uma radiação Hawking acústica, na forma 
de um banho térmico de fônons, que é caracterizada pela temperatura acima. 

É importante observar que, nossa métrica acústica depende da densidade do 
fluido, da velocidade do fluxo e da velocidade local do som. Ela é governada pelas 
equações de movimento do fluido somente, não utilizamos em momento algum, as 
equações de Einstein da relatividade geral. Ou seja, nosso modelo só pode ser usado 
para uma analogia com a cinemática do sistema gravitacional. Sendo os aspectos 
dinâmicos não inclusos neste modelo, ficamos impossibilitados de reproduzir a 
termodinâmica de um buraco negro. Este resultado não é restrito ao nosso modelo 
acústico, vários outros modelos que se propõem a imitar a gravidade não estão 
sujeitos ã dinâmica do relatividade geral, veja por exemplo \27]. 
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No entanto, vamos ver no próximo capítulo que, para um buraco negro 
bidimensional, podemos sim obter uma definição de massa para o buraco negro 
acústico e assim produzir uma analogia às leis da termodinâmica do sistema 
gravitacional. 
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Capítulo 5 

Analogia à dinâmica de um buraco 
negro 

Nas sessões que seguem este capítulo iremos demonstrar que, a partir de um 
buraco negro bidimensional, podemos obter uma analogia direta entre a dinâmica 
do sistema gravitacional e a dinâmica do fluido. Veremos que reduzindo um buraco 
negro 4D esfericEimente simétrico a uma descrição 2D, sua dinâmica equivale a uma 
dinâmica de fluido com um certo vínculo, que nos permitirá definir a primeira lei da 
termodinâmica para um buraco negro acústico. 

5. 1 Métrica acústica para um buraco negro 2D 

Um sistema gravitacional quadri-dimensional esfericamente simétrico, tal como o 
espaço-tempo de Schwarzschild, pode ser representado por um sistema bidimensional 
com um campo dilatônico caracterizado por um potencial V, cuja ação é dada por 
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S=^j(fxV^{^R + X^V{4>)), (5.1) 

onde (p é um campo escalar, o próprio dílaton, e A é um parâmetro com unidade de 
[comprimento]^^. Aqui estamos usando c = h = kB = l. Esta ação admite soluções 
tipo Schwarzschild que toma a forma 



ds' = -( J(0) - — J dti + (^J(</>) - — J dH, (5.2) 



onde (/) = Ar, M é a massa do buraco negro e J = J Vdcj). O horizonte de eventos ocorre 
quando J{4)h) = 2M/A e assim (ph = Ar/^. Agora precisamos encontrar uma métrica 
acústica análoga a equação j5.2| . 

Do capítulo anterior temos que, a métrica acústica para um fluido ideal de fluxo 
unidimensional, é escrita como 



ds^ = ^ [-{cl - vl)df - 2vodxdt + dx^] , (5.3) 



e a velocidade local do som é dada por = dP/d^. Usamos para identificar a 
densidade usual, pois iremos trabalhar com uma densidade adimensional po = A^^põ- 
Da mesma maneira como foi feito anteriormente, devemos fazer uma transformação 
de Painlevè-GuUstrand na métrica da equação |5.2 para fins de comparação com 



a métrica acústica. Usando a equação (4.20| , que define um perfil de velocidade 
para o qual o sistema acústico é conforme ao sistema gravitacional, reescrevemos a 
transformação j4. 18) como sendo 



dts = dt+ „ dx. (5.4) 
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Substituindo isso na métrica j5.2) , teremos 



ds' 



J 



2M 



J r- ) — ^dtdx 



J 



2M\ vá 



2M 



, 1 ,^^,dx^ + U - ^ 1 dr\ (5.5) 



-1 



Comparando a componente tempo-tempo dessa métrica com a métrica acústica usual 
escolhemos que 

(5.6) 



2M «0.2 2^ 

J = — + —{c,-Vo), 



A c. 



e então obtemos 



ds 



2 Po 



(c^ — VQ)dt'^ — 2vodtdx ^ — ^^—^dx 



o 



+ 



(5.7) 



Por fim, façamos uma transformação na coordenada radial do sistema gravitacional 
para a coordenada cartesiana do sistema acústico, 



dr = podx, 



(5.8) 



e então 



ds' 



Po 



{cl — VQ)dt' — 2vodtdx 



{cl - vl) 



-dx^ + 



{cl - vl) 



dx' 



ds^ = — [-{cl - vl)dt' - 2vodtdx + dx^] . 

Cg 



(5.9) 



Assim sendo, vemos que, com as transformações dadas por |5.4 , |5.6 e |5.8 



, a 



correspondência entre os dois modelos é exata. Isso difere do modelo análogo para 
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um buraco negro 4D, como visto no capítulo anterior, onde a equivalência entre eles 
só é válida perante transformações conforme da métrica. 

No entanto, este modelo é limitado a imitar somente a cinemática de um buraco 
negro 2D. Para que possamos estender a analogia até um nível dinâmico é preciso 
comparar a dinâmica gravitacional com a dinâmica do fluido, que determina po^^^o e 

Cs. 

5.2 Analogia à dinâmica gravitacional 

A principal diferença entre os dois sistemas físicos tratados nesta dissertaçáo 
é que o sistema gravitacional deve ser covariante sobre transformações gerais da 
métrica. A relatividade geral é uma teoria independente do "background", ou seja, 
não é necessário a existência de uma métrica específica para deduzirmos as equações 
de campo provenientes da ação |29[ , |30| . Já no sistema acústico a métrica surge de 
perturbações nas equações de movimento, sendo então já previamente fixada. 

Portanto, para construir uma dinâmica acústica análoga, devemos fixar uma 
métrica para a descrição do nosso sistema gravitacional dilatõnico 2D. Porém, a 
gravidade dilatõnica 2D é puramente topológica, então não temos uma propagação 
física dos graus de liberdade gravitacionais. Isto deixa em aberto a possibilidade de 
que, uma vez que tenhamos fixado a métrica, a dinâmica dos graus de liberdade 
gravitacional possa ser imitada pela dinâmica do fluido. 
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As equações de campo para o modelo dado em |5. 1) são 



fí = -A^, (5.10) 
V^V,0 - = 0. (5.11) 



Nós fixamos a métrica escolhendo o calibre de Schwarzschild 



ds' = -X{r)dti + X{r)-'dr^, (5. 12) 

que é estática, justamente como queremos para que seja compatível com o modelo 
acústico já apresentado aqui. Utilizando as equações \5. 1 1 e (5. 12 obtemos que 



^=A, (/) = Ar, (5.13) 
dr 



^ = AF, (5.14) 

dr 



de onde pode-se verificar que o buraco negro (5.2 é solução dessa última expressão 



Isso nos mostra que, quando fixamos o calibre, o sistema gravitacional fica com 
somente um grau de liberdade, a função X{r) que parametriza a métrica, e uma 
única equação. O campo escalar (p tem como vínculo ser dtretamente proporcional a 
r, portanto ele representa apenas uma coordenada tipo espaço. 

Agora precisamos obter as equações da dinâmica dos fluidos, que são as equações 
da continuidade e de Euler, respectivamente representadas pelas equações (4.1 e 



(4.2) . Mas faremos um tratamento unidimensional. Assim, as equações anteriores se 
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resumem a 



dvQ dp _ dip ^ ,^ , ^, 

ax dx ax 



po{x)vo{x)A{x) = D, (5.16) 

onde A denota a área da secção transversal por onde este fluido está passando e 
D é uma constante. Vamos considerar dois casos diferentes para solucionar estas 
equações: 

a - O potencial externo i) é não homogéneo, sendo dependente da posição, e ^4 é uma 
constante; 

b - O potencial externo é nulo e Aé uma função dependente da posição x. 
Caso - a 



Vamos reescrever as equações |5.15) e |5. 16) a partir da mudança da coordenada 
pela coordenada radial do sistema gravitacional tal como foi feito antes usando (5.8 



Assim, usando também a relação entre e p^, e considerando a condição baro trópica 
que implica que a pressão só dependa da densidade do fluido, temos que 



dvo dp dpo dtp 1 71 

ar apo ar ar 



Multiplicando tudo por 2/cs e usando a definição de velocidade do som onde = 
dp/dpo, temos que 

2Po^d^^2c.^ + 2^? = 0. (5.18) 

dr dr dr 
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Definindo então as novas variáveis 



Y = pocs, (5.19) 

F = In ' 

■ PO 



obtemos 



dZ ^ 1 dpo 2 ^,2n 2 ^°^° 

dr dr " ^ c, dr ' 



dy dpo 
dr * dr ' 



dF 1 dyOo 

dr Po dr ' 



Co 



para reescrever a equação de Euler como 



Na equação (5.16) escrevemos que 



dr 



dr 



(5.20) 



ÍE=2í5--zíí+2.-* (5.211 

dr dr dr dr 



(/9ofo) = O, (5.22) 
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mas como 



logo 



Pfívl = pocl - CsX poVo = VYiy - Z), (5.23) 



^ [Y{Y - Z)] = 0. (5.24) 



Com as equações j5.24) e |5.21) representando a dinâmica do fluido vemos que, 
como a analogia entre as métricas é feita tal que Z = X (componentes temporal de 
cada métrica) , a equivalência entre a dinâmica gravitacional e as equações de Euler e 
da continuidade pode ser estabelecida impondo um vínculo de acordo com a equação 
(5.14) . Ou seja 



dY dF F dé , 

!— - Z— + 2e-^ = XV. (5.25) 
ar ar dr 



Isso mostra que, se o potencial externo for dado, a dinâmica do fluido é 
caracterizada por três funções desconhecidas, po,vo,Cs e duas equações. Temos então 
a liberdade de estabelecer um vínculo como uma terceira equação para o sistema. 
Na verdade a equação (5.25 fornece uma equação de estado para o fluido para um 



dado ^/J, e assim esta equação de estado pode ser vista como uma consequência 
da dinâmica gravitacional. Por outro lado, se tivermos uma equação de estado 
previamente definida para o fluido, este vínculo definirá o potencial externo. 

A massa para o buraco negro dilatõnico em 2D fST 



M = ^ y V{cP)dcP - (Vcpf^ , (5.26) 

permanece constante quando calculada no horizonte de eventos. Seguindo esta 
formulação, podemos obter uma expressão equivalente a massa do buraco negro de 
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Schwarzschild através das equações j5.24 e j5.21| que também seja uma constante 
independente da coordenada espacial x. Definimos então a função F, tal que 



dT 
dr 



dY 
) 

dr 



dr 



dr 



dZ' 

dr 



(5.27) 



Portanto, escrevemos que 



T = a 



-ZdF + 2e-^ dij ] + 2Y - Z 



(5.28) 



onde a e ro são constantes arbitrárias. Porém, podemos determinar um valor para a 
constante a tal que T = M e usando Z = J — 2M/X, assim a = A/2. 

Com a função T fazendo o papel da massa do buraco negro acústico, se obtivermos 
loma expressão para a temperatura Hawking equivalente, podemos formular as leis 
da termodinâmica para o modelo análogo. Assim, usando a equação j3.23) para a 
métrica acústica, obtemos a gravidade superficial acústica 



1 dZ 

2 dr 



r=rh 



(5.29) 



e da definição da temperatura Hawking, T = k/2tt, temos que 



Ta 



1 dZ _ 1 
47r dr Att 



Cs=VO 



(5.30) 



que é exatamente o que obtemos no capítulo anterior. Podemos escrevê-la também 
em função das variáveis novas. 



1 



dV dF 



i2— -Z— + 2e- 
47r/9o \ dx dx 



dip 
dx 



(5.31) 
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Para a massa Ma do buraco negro acústico, tomamos a equação j5.28 calculada 
no horizonte de eventos com a = X/2. Então, 



í dF p dij)\ , : 



(5.32) 



onde usamos que Z{xh) = 0. Com isso podemos escrever que 



dMa = Ta{2TTXpodx). (5.33) 



Definindo a função 



Sa = 2ttX í podx, (5.34) 

J oo 



obtemos finalmente o equivalente a primeira lei da termodinâmica, 

dMa = TadSa, (5.35) 

onde Sa é a entropia do buraco negro acústico. É importante notar que, se tivéssemos 
utilizado diretamente a definição de entropia do sistema gravitacional contido nas 
equações (3.56 e |3.57) , teríamos obtido o mesmo resultado para a entropia do 



sistema acústico. Isso equivale dizer que Sa é equivalente a S até mesmo em sua 
definição. 

Caso - b 

A seguir mostramos que os resultados são exatamente os mesmos quando 
consideramos o potencial externo como sendo nulo e adotamos a área da seção 
transversal por onde o fluido escorre como imia fiinção da posição. Neste caso, as 
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equações da dinâmica do fluido escritas em função das variáveis (5. 19 , são 



dZ _^dY ^dF 
dr dr dr ' 



d^ 
dr 



[Y{Y - X)A^] = 0. 



(5.36) 



Da mesma forma como fizemos no caso anterior, introduzimos o vínculo 



ar dr 



(5.37) 



e assim as equações j5.36| podem ser feitas equivalentes com a dinâmica gravitacional 
da equação |5.14) . Nesse caso, o vínculo acima também representa uma equação de 
estado para o fluido se A{x) for dado. Por outro lado, se a equação de estado for 
conhecida, o vínculo representará uma expressão para A{x). 



A função constante que pode ser interpretada como a massa do buraco negro 
acústico é 



V = a 



/ {-ZdF) + 2Y - Z 



E desta equação obtemos a massa Ma 



(5.38) 



Mn 



dx 



XQ 



(5.39) 



e a temperatura Hawking acústica 



2^^ - 

dx dx 



(5.40) 



Assim a primeira lei da termodinâmica é respeitada com a entropia do sistema 
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acústico dada novamente pela equação j5.34) . 
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Considerações Finais 



Este trabalho apresentou uma revisão teórica sobre um modelo acústico análogo 
ao buraco negro de Schwarzschild. Explicitamos a maneira como esta analogia pode 
ser estabelecida quando perturbamos as equações da dinâmica dos fluidos e com isso 
conseguimos descrever a propagação das ondas sonoras no ferramental matemático 
da relatividade geral de Einstein. Essas ondas sonoras, que foram tomadas como 
sendo pertubações no campo de velocidade de propagação do fluido, descrevem 
geodésicas nulas em um espaço-tempo efetivo governado por uma métrica dependente 
da velocidade local do som, da densidade e da velocidade do fluido: 



„2 7.2 



dr + 5ij{dx^ - Vodt){dx^ - v^dt) 



(6.1) 



A métrica acústica obtida tem termos cruzados do tipo tempo-espaço e, para 
efeito de comparação, devemos escrever a métrica de Schwarzschild na forma 
de Painlevè-Gullstrand (equação |4.19) . Notamos que a equivalência entre essas 
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duas métricas é conforme, o que a princípio não apresenta problemas quanto 
ao uso dessa analogia para o estudo da radiação Hawking visto que a gravidade 
superficial e, consequentemente, a temperatura Hawking são grandezas invariantes 
sob transformações conformes da métrica. 

Para obter um equivalente acústico da gravidade superficial de um buraco negro 
procedemos da mesma maneira que foi apresentado no capítulo 2. Através da 
definição de horizontes de líilling obtém-se uma relação entre a gravidade superficial 
e a derivada dos elementos da métrica. O resultado obtido é exatamente o mesmo que 
aqueles apresentados em )18[|21), mas que foram calculados de uma forma distinta. 



1 d 



1 d f 2 2^ 



(6.2) 



Sobre o fato de obtermos uma temperatura Hawking para o buraco negro 
acústico, equação |4.28 , conclui-se que a radiação Hawking é um fenómeno 



puramente cinemático jã que não utilizamos as equações de campo de Einstein 
para descrever o nosso modelo |33) . Isto é, para uma geometria Lorentziana, basta 
que haja lom horizonte de eventos que a radiação Hawking surgirá sem que seja 
necessário previamente impor uma dinâmica para o sistema. No entanto, o fato do 
modelo acústico ter uma métrica conforme e não exatamente análoga ã métrica de 
Schwarzschild, gera complicações se tentarmos imitar a dinâmica da relatividade 
geral com este modelo. Mais especificamente, uma analogia ã entropia de um 
buraco negro seria bastante complicado pois teríamos de fazer modificações na ação 
da mecânica dos fluídos, já que esta entropia está diretamente relacionada com o 
teorema da área do horizonte de eventos de um buraco negro, explicitado na equação 
{3:56) . 
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Mesmo assim, essa dificuldade pode ser vencida se trabalharmos em cima de uma 
analogia a um buraco negro bidimensional. Este sistema 2D é obtido utilizando 
a descrição do buraco negro de Schwarzschild através de um campo escalar, o 
dílaton. Com isso, a equivalência entre as métricas ocorre de forma exata e o sistema 
acústico é somente bidimensional também, permitindo estender a analogia até o nível 
dinâmico. 

Neste ponto fica claro uma diferença crucial entre os modelos gravitacionais e 
actistico. A relatividade geral é uma teoria independente da métrica, as equações de 
campo são derivadas de uma determinada ação não sendo necessário a definição de 
imia métrica para isso. Já no modelo actistico, a métrica que governa a propagação 
das ondas sonoras é obtida na pertubação das equações de Euler e da continuidade, 
ou seja, uma consequência que surge quando descrevemos a propagação do som 
usando o ferramental matemático da relatividade geral. Isso implica que devemos 
previamente impor um vínculo no sistema gravitacional fixando uma métrica para 
ele, neste caso uma métrica no calibre de Schwarzschild. 

ds^ = -X{r)dtl + X(r)-^dr^ (6.3) 

Feito isso podemos definir uma função equivalente à massa do buraco negro 
gravitacional e, a partir dessa função, obter uma. entropia para buraco negro acústico 
assim como a primeira lei da termodinâmica do buraco negro acústico. No entanto 
este modelo também apresenta suas limitações devido as aproximações utilizadas 
sobre ele. Essa analogia ã dinâmica de um sistema gravitacional só é válida para 
soluções estáticas, o que nos impede de utilizá-las em várias situações interessantes 
onde o caso não estático é essencial, por exemplo, uma investigação da evaporação 
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de um buraco negro através do cálculo das reações de fundo, fenómeno que nem 
mesmo Hawklng levou em consideração quando realizou seus cálculos. Também há 
o fato de náo considerarmos um modelo 2D com massa. A ausência de fontes para o 
campo gravitacional representa um enorme simplificaçáo no modelo pois as equações 
dinâmicas ficam muito mais simples. Mas, sabemos que é nescessário a presença de 
fontes para descrevermos uma situaçáo realística, porém, encontrar uma grandeza 
no sistema acústico similar ao tensor momento-energia do campo gravitacional pode 
se tornar algo muito complicado de se fazer. 

Como uma sequência desta dissertação, pretende-se aplicar simulações numéricas 
ao estudo da física dos buracos negros acústicos ( |34| , (35), (36)) náo restringindo 
a análise a uma regiáo muito próxima do horizonte sónico, isto é, a regiáo onde 
o movimento do fluido muda de subsónico para supersônico, em conexáo com 
modelos de Física da Matéria Condensada. Neste sentido, tal procedimento aplicado a 
correntes de fluido (inclusive correntes de fluido em bocais de Lavai) pode averiguar a 
formação de buracos negros acústicos e estimar a contrapartida clássica da radiação 
Hawking, ou seja, o espectro de potência da perturbação do potencial velocidade do 
fluido. 
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